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We investigate the nature of the new cross term for Gaussian 
parameterizations of Bose-Einstein correlations of identical 
particles emitted from purely chaotic hadron sources formed 
by relativistic heavy ion collisions. We find that this addi- 
tional parameter in the so-called Bertsch parameterization 
can be expressed in terms of a linear "out-longitudinal" cor- 
relation coefficient for emission of bosons and two already 
known "radius" parameters, Ri and R . The linear corre- 
lation coefficient is of kinematical nature and can be used 
to determine the widths of longitudinal momentum distribu- 
tions. 

In the ongoing search for the quark-gluon plasma (QGP) 
one is especially interested in the volumes and the life- 
times of the hot and dense zones of nuclear matter, i.e., 
the fireballs that are generated in relativistic heavy-ion 
collisions. A tool to measure space-time extensions || of 
the fireball is given through the Hanbury-Brown/Twiss 
(HBT) interferometry g], which provides the measure- 
ment of Bose-Einstein correlations (BEC) Q. 
BEC functions of identical bosons are quantum- 
statistical observables which contain information about 
ten quantities which characterize a hadron source: 
lifetimes, longitudinal and transverse extensions of the 
chaotic and the coherent source, temporal and spatial 
(longitudinal and transverse) coherence lengths and the 
chaoticity. In the case of pion interferometry it has been 
argued that a possible coherent subcomponent of the 
hadron source is almost unobservable flq], because reso- 
nance decays contribute strongly to BEC. Many param- 
eterizations of experimental BEC functions are based on 
a neglect of a coherent component thereby resulting in 
the description of the hadron-emitting source with only 
three quantities, (e.g., lifetime, longitudinal and trans- 
verse radii of the purely chaotic hadron source). 
For heavy-ion collisions a very popular parameteriza- 
tion of purely chaotic BEC functions is provided through 
the so-called Bertsch parameterization |6y wherein a two- 
particle BEC function is parameterized in terms of the 
three momentum-dependent "radius" parameters Ri(K), 
R S (K), and R (K). Here, K = ^(ki + k 2 ) is the aver- 



age momentum of two identical bosons which are emitted 
with their individual momenta, fei and k 2 , respectively. 
In the HBT Cartesian coordinate system the "longitudi- 
nal" or z direction (subscript I) is parallel to the beam, 
the "out" or x direction (subscript o) is parallel to the 
component of K which is perpendicular to the beam, 
and the "side" or y direction (subscript s) is the remain- 
ing direction. In ref. |7]] it was shown that even in the 
case of a purely chaotic hadron source, three quantities 
are not adequate to fully parameterize the system under 
consideration. E.g., a cylindrically symmetric expanding 
fireball has in general to be described with an additional 
"cross-term" radius[], R i(K). 

The generalized (Gaussian) Bertsch parameterization of 
a purely chaotic two-particle Bose-Einstein correlation 
function is given through J7| 



C 2 (ki,k 2 



1 + \(K) ■ exp[- qfRf(K) - q 2 R 2 (K) 
-q 2 R 2 s (K) -2q oqi R 2 ol (K)}. 



(1) 



In eq. (|l|) the qi (i = I, s,o) refer to the components of 
the momentum difference q — k\ — k 2 , and X(K) is the 
momentum dependent incoherence factor which accounts 
for reductions of the BEC due to long-lived resonances 
|| and averaging due to phase-space, respectively. 
In a physical interpretation one is troubled by the fact 
that one and the same two-particle BEC function is de- 
scribed in case (a) of the generalized Bertsch parameter- 
ization through four "radius" parameters (Ri, R s , R , 
Roi), and (b) in case of the Yano-Koonin-Podgoretskii 
(YKP) parameterization |j] through three "radius" pa- 
rameters (i?0j Rh R±) an d one kinematic quantity, the 
Yano-Koonin velocity (v). It is the purpose of this paper 
to take a closer look at the recently reported new cross 
term, R i, and to investigate its nature. 
Before we begin our analysis, we recall that the four 
"radius" parameters of the generalized Bertsch param- 
eterization can be expressed in terms of variances of 
the (boosted) components of the space-time points x v = 
(t, x, y, z) G S (E is the freeze-out hypersurface) and co- 
variances between the components of x v (cf. ref. B) 
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1 Experimentally, a "cross-term" has been observed, e.g.. 
the NA35/NA49 Collaborations pi. 



by 



R 2 S {K) « a 2 y , 
Rf(K) 



'z-p«t 



R 2 o(K) 

R 2 ol (K) * o 



'x- fixt ' 



(2) 



x — [3±t, z — j3\\t • 



In eg. (ph w e use $ = KiJEk (i =||,-L), where Ek = 
V m 2 + K 2 . 

Let us now consider a quite general mathematical repre- 
sentation of a two-dimensional Gaussian, f(gi,q ), in the 
two variables qi and g with fixed parameters Ri, R 0l and 
A 

f(q h q ) = 1 + A • cxp[- q 2 RJ - q 2 R 2 

+ 2 q qiR R[ p i] , (3) 

where p i is the "linear correlation coefficient" which 
gives the strength of linear correlation between the two 
quantities qiRi and q R . A comparison between eq. (|l|) 
and eq. (f|) yields for the "cross-term" (while ignoring 
the "side-term" without loss of generality) 



Rli(K) 



p ol (K) ■ R (K)Ri(K), 



(4) 



or in terms of the variances and covariances for the 
(boosted) "longitudinal" and "out" components of x v 



Poi(K) = - 



RliiK) 



Ro{K)Ri{K) 



&x-0xt, z-f)\\t 

2 2 

a x-0j_t a z-f3 u t 



(5) 



with \p i(K)\ < 1 Vl, representing the Cauchy- 
Schwarz inequality relation for p \. 

Thus the "cross-term" radius parameter, R i in the 
Bertsch parameterization can be expressed in terms of a 
linear "out-longitudinal" correlation coefficient, p i, for 
emission of bosons and the two already known "radius" 
parameters, Ri and R . 

In order to illustrate the nature of the cross-term, in Fig.l 
we plot the pion rapidity distribution (top) and the scat- 
ter plot qi vs. q (bottom). These plots are the results 
of the microscopic transport model RQMD (l0| calcula- 
tions for 200 AGeV S+S central collisions. All plots are 
generated in the nucleon-nucleon center of mass system; 
the rapidity windows are labeled in the figure. 
In the mid-rapidity window no correlation between qi 
and q is seen while sizable correlations are found in 
other rapidity windows. In addition the strength of the 
correlation depends on the rapidity for which one con- 
structs the scatter plotQ Qualitatively, the strength of 



2 Note that in case of a purely Gaussian distribution a non- 
vanishing p i is represented by a tilted ellipsoid. However, 
inspecting Fig. 1 we found a non-ellipsoid shape in the scatter 
plots. This feature indicates non-Gaussian distributions in 
real collisions. This will not alter our main conclusion in the 
text. 



the correlation is consistent with the NA35/NA49 result 
|g|. No correlation is seen between qi and q s or q s and 
q in this model calculation. As discussed later in more 
detail, this longitudinal momentum dependence implies 
that the cross-term is kinematic in nature. However, as 
can be seen in the figure this dependence could be used 
to determine the rapidity width of the particles under 
study. 

For the sake of further illustration we consider results of 
a hydrodynamical analysis 0@ of 160 AGeV Pb + Pb 
(central) collisions using the HYLANDER code pi. In 
the following we discuss the fireball emission of (identi- 
cal) K~ K~ pairs because in ref. |jj it has been argued 
that for K~K~ BEC functions their numerical calcu- 
lation and subsequent fit can be replaced by the direct 
evaluation of eqs. (Q) in order to determine the four "ra- 
dius" parameters of the generalized Bertsch parameteri- 
zation. We perform our analysis for average, transverse 
kaon pair momenta K± = 1.5 GeV/c (an analysis for dif- 
ferent average transverse kaon pair momenta, K± =/= 0, 
leads qualitatively to the same results). 
We now introduce the strengths of linear correlation 
among t, x, and z 



Pt,x 



&t,x 



yM 



Pt,z 



<?l. 



yff 



Px,z 



\A 



(6) 



From eqs. 
and R i 



and (o) we are now able to rewrite Ri, R , 



R 2 (K) » ul + Pla 2 - 2 (3 ± a t a x p t>x , (7) 

Rf{K) w o\ + 0\a\ - 2 f3 l{ a t a z p M , 
Rli(K) w a x a z p XtZ + )3\\P±a 2 

- /3\\(Ttcr x Pt,x - (3_Lcr t <T z p t , z . 

In Fig. 2 we show the linear correlation coefficients due 
to eqs. (o) and (0) for emission of negative kaons as 
a function of the kaon-pair rapidity, y^, at K± = 1.5 
GeV/c. The results have been obtained for Pb + Pb at 
160 AGeV with a rapidity, t/a = 0.90, at the edge of 
the initial central fireball (cf. ref. |12] for more details). 
Fig. 2 is divided into three ranges of the kaon-pair ra- 
pidity, hk < — l-4, \uk\ < 1-4 ,and 1.4 < i/k, which 
correspond to the freeze-out hypersurface emission zones 
SWR, NTR, SWR, respectively. SWR denotes the "sim- 
ple wave region" which is given by the original Riemann 
rarefaction waves, and NTR denotes the "non-trivial re- 
gion", where two rarefaction waves overlap (cf. [Q and 
also Fig. 3). I.e., kaon pairs with \yx\ < 1.4 are mostly 
emitted from the NTR, while kaon pairs with 1.4 < \yx\ 
are mostly emitted from the SWR. 

Fig. 2(a) shows that there is an almost perfect linear 
correlation between the emission times t and the trans- 
verse emission coordinates x (the longitudinal emission 



coordinates z) in the NTR (SWR) of the freeze-out hy- 
persurface. We have \pt,x\ ~ 0.85 (dotted line) in the 
NTR, and \p t . z \ « 1.0 (dashed line) in the SWR. For 
\i/k\ — ► oo we get \pt,x\ — ► 0, and because of symmetry 
reasons we have pt,z{yK = 0) = 0. The linear correla- 
tion coefficient p x _ z (solid line) reaches its extrema at the 
transition zones of NTR and SWR. 

Fig. 2(b) shows that by coincidence p x%z w p tx ■ p tz , but 
this relation is in general not true. Fig. 2(c) shows that 
there is a direct proportionality between the (internal) 
linear correlation coefficient, p x>z (solid line), and the 
measurable p i (dashed line): 



Poi(K) ex p x 



(8) 



Because of this proportionality p i reflects the strength 
of linear correlations between transverse and longitudi- 
nal emission coordinates of the freeze-out hypersurfacc. 
A negative value for p \ at positive average kaon-pair ra- 
pidities indicates that for larger longitudinal emission po- 
sitions, the transverse emission positions become smaller. 
Conversely, a positive value for p i at negative average 
kaon-pair rapidities indicates that for larger longitudinal 
emission positions (in the negative direction) the trans- 
verse emission positions also become smaller. Hence, p i 
enables us to measure finite-size effects of fireballs in ra- 
pidity space. 

In Fig. 3 we show (a) the rapidity spectra, dN/dy, (b) 
linear correlation coefficients, p i, and (c) freeze-out hy- 
persurfaces S(r = 0) for emission of negative kaons for 
Pb + Pb at 160 AGeV. The results have been obtained 
with different rapidities, j/a, at the edge of the initial 
central fireballs. The parameter j/a controls the strength 
of the initial rapidity field of the initial fireball in the hy- 
drodynamical model. A larger value for j/a results in a 
larger initial longitudinal rapidity field of the fluid, and 
subsequently at freeze-out it leads to rapidity spectra of 
larger widths (c/. Fig. 3(a)). Fig. 3(c) shows the corre- 
sponding freeze-out hypersurfaces £(r = 0). Consistent 
with expectation, a larger value for y\ results in more 
longitudinally extended freeze-out hypersurfaces. Fur- 
thermore, we observe that with increasing j/a, the NTR 
of the freeze-out hypersurfaces becomes more and more 
hyperbola shaped (like the scaling solution for Bjorken 
initial conditions fl5[| ). 

We stress, as shown with Fig. 3(b), that the Full Widths 
at Half Maximum of the rapidity spectra, dN/dy, of neg- 
ative kaons nearly coincide with the positions of extrema 
of the corresponding linear "out-longitudinal" correlation 
coefficients, p i{K). We also observe that a more hy- 
perbola shaped freeze-out hypersurface results in smaller 
absolute values for p i. 

After introducing the linear correlation coefficient, p i, 
the difficulties of interpretation encountered above no 
longer persist, if one uses in case of the generalized 
Bertsch parameterization the three "radius" parameters, 



Ri, R s , and R , and the kinematical quantity p i. In gen- 
eral, if one analyzes a heavy ion collision which has no 
azimuthal symmetry, further (measurable) linear correla- 
tion coefficients come into existence. In case of the gen- 
eralized (Gaussian) Bertsch parameterization, one would 
use for a fit of a two-particle BEC function (i,j — 1,8,6) 



C(K,q) = l + \{K) 



(9) 



x exp 



J2 QiRt(K) QiRiiX) ~Pij{K) 



with 



Pij(K) 




1 



d 2 C(K,q) 



2Ri(K)Rj(K) dqtdqj 

®Xi-f3it, Xj-fjjt 



9 = 



2 2 

G Xi-(3it (7 x j -(3jt 



(10) 

(ii) 

(12) 



and \pij(K)\ < 1 \/K. In cq. (JI4), we have used the 
notation xi = z, x = x, x s = y, Pi = f3\\, (3 = f3±, and 
& = 0. 

To summarize, we have investigated the cross term for 
Gaussian parameterizations of Bose-Einstein correlation 
functions of purely chaotic hadron sources. We find that 
at relativistic energies this additional parameter in the 
so-called Bertsch parameterization can be expressed in 
terms of a linear "out-longitudinal" correlation coeffi- 
cient, p i(K), for emission of bosons and the two already 
known "radius" parameters, Ri and R . The correlation 
coefficient provides an additional tool for determining the 
particle rapidity distributions. 
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FIG. 2. Linear correlation coefficients calculated from HY- 
LANDER, for emission of negative kaons for Pb + Pb at 160 
AGeV as a function of the kaon-pair rapidity, ypc, at K± = 
1.5 GeV/c. (a) p t , x (dotted line), p t , z (dashed line), p x . z 
(solid line); (b) the product p t , x ■ pt.z (dashed-dotted line), 
px,z (solid line); (c) the observable p oi (dashed line), and p x , z 
(solid line). 



FIG. 1. Rapidity spectrum of negative pions in arbitrary 
units for S + S at 200 AGeV from RQMD (version 2.1). Pion 
pairs have been formed in different rapidity bins and their cor- 
responding correlations between q and qi are shown as scatter 
plots. The different colors in these plots indicate different rel- 
ative strengths of the correlation functions (seven equidistant 
regions). The correlation functions take their maximum val- 
ues (near to 2) in the center of each scatter plot. The large 
red areas indicate a minimum correlation strength near to 1. 



FIG. 3. (a) Rapidity spectra, dN/dy, (b) linear correlation 
coefficients, p i, at K± = 1.5 GeV/c and (c) freeze-out hyper- 
surfaces E(r = 0) for emission of negative kaons for Pb + Pb 
at 160 AGeV. The four line styles correspond to four different 
rapidities, yA, at the edge of the initial central fireballs. 



